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ORBIT DECAY ANALYSIS OF 
STS UPPER STAGE BOOSTERS 

by 

Otis F. Graf, Jr. and Alan C. Mueller 


1.0 BACKGROUND AND INTRODUCTION 


During the next decade, the Space Transportation System 
(STS) will launch many missions to geosynchronous orbit. 
Expendable upper stages will be used to boost the payloads 
out of near earth orbit and into the geosynchronous transfer 
orbit. After the mission is completed, these boosters re- 
main in the transfer orbit. It is of some concern that over 
the years, these upper stage boosters could pose a hazard to 
future space missions. Thus, the purpose of this task was 
to study the lifetimes of the geosynchronous transfer orbits 
and to analyze the orbit decay. 

This task (Exhibit "D" , Contract NAS9-15444) is a follow- 

f 

on to a previous study under Contract NAS9-15444. That study 
showed that the geosynchronous transfer orbits can experience 
large variation in perigee altitude. It was also shown that 
the orbitswill rapidly decay for some initial conditions. 
However, the previous studies were based on a trajectory pro- 
gram (DSTROB) that did not include the effects of atmospheric 
drag. Thus,' actual orbit decay studies were not carried out 
in the previous study. 


This study is reported in ACM Technical Report TR-115 
"A Study of the Lifetimes of Geosynchronous Transfer 
Orbits" . 


The next section contains a summary of the work that was 
carried out under Exhibit "D". Appendix A contains the slides 
that made up the Final Presentation of the work to the Johnson 
Space Center. These slides include details on the launch 
window studies. Appendix B is a copy of a paper that was pre- 
sented at the 1979 AAS/AIAA Astrodynamics Specialists Conference. 


2 . 0 SUMMARY 


»J* 

Atmospheric drag effects were added to the DSTROB 
trajectory program. The perturbation models are described 
in Appendix B. The program was optimized for fast execution 
in a demand (interactive) environment. A twenty year time 
history of a geosynchronous transfer orbit can be computed 
in about 5 seconds . The accuracy of the DSTROB algorithm 
was checked against precision numerical integration methods. 

Orbit lifetimes of the transfer orbit were studied for 
several planned STS geosynchronous missions . Included in 
the study were the SUSS-A and SUSS-D boosters on STS-G, 
the IUS booster on STS-7, and the SUSS-A booster on 
STS-8. Time history plots of perigee altitude were gener- 
ated, using the DSTROB program. These plots are included 
in Appendices A and B. 

Because of the fast execution time of DSTROB , it can 
be effectively used for parametric studies and scans. A 
driver routine for DSTROB was developed, that allowed orbit 
lifetime plots to be generated as a function of the initial 
longitude of the ascending node. This driver routine will 
generate trajectories on 10° increments of ascending node, 
and store the lifetime values on a mass storage file. The 
data can then be submitted to a plot routine for automatic 
output of figures. It was found that the initial ascending- 
node and epoch are important parameters in determining orbit 
lifetime . 

The study also included development of launch windows, 
such that the transfer orbit will decay within two years. 
Launch windows were given as a function of either initial 


This is prototype program that was developed under 
Exhibits "C" and "D" of Contract NAS9-15444. DSTROB 
resides on the JSC Univac 1108/8 system. 




ascending node or launch time. Example launch windows are 
included in Appendix A. These launch window figures were 
drawn by hand, using data that was generated by DSTROB , 
and a calculator program to compute launch time. However, 
the process could be automated so that a computer program 
would produce a hard-copy launch window figure. 
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OVERVIEW OF THE DSTROB ALGORITHM 

DISCUSSION OF THE ATMOSPHERIC DRAG 
MODEL AND COMPARISON RESULTS 

ORBIT LIFETIMES OF STS~6 AND STS-7 
PAYLOADS 

SOME ANALYTICAL RESULTS 

LAUNCH WINDOWS DERIVED FROM ORBIT 
DECAY STUDIES 

SUMMARY AND CONCLUSIONS 




ORBIT GEOMETRY 


SUMMARY OF THE STROBOSCOPIC METHOD 


DEFINE 


WHERE a 


n 


Aa n = a„ - a n 
n n + i n 

IS THE VECTOR OF ELEMENTS AT THE 


11 “TH PERIGEE. 


Alf n IS GIVEN BY AN ANALYTICAL EXPRESSION: 


SUMMARY OF THE STROBOSCOPIC METHOD 
(concluded) 




AERODYNAMIC MODEL IN DSTROB 

• CONSTANT BALLISTIC NUMBER VALUES USED IN CURRENT 
SIMULATION: 

C (1 * 2.2 


I US 

MASS 65 2200 lbs. (098 kg.) 
AVERAGE AREA - 84 it 2 (7,« „i 2 ) 

SUSS-A 

MASS * 700 lbs, ( 316 kg,) 
AVERAGE AREA - 42.2 ft 2 (3.02 in 2 ) 

SUSS-D 

MASS - 515 lbs, (234 kg.) 
AVERAGE AREA - 27.4 ft 2 (2,55 in 2 ) 






m 
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AERODYNAMIC MODEL IN DSTROB 

(concluded) 


• DENSITY MODEL 

P = Pp exp 

P i DENSITY AT PERIGEE 

r 

h : PERIGEE ALTITUDE 

r 

h g s SCALE HEIGHT 
h : CURRENT ALTITUDE 

COMMENTS l 

O) p n AND h e ARE COMPUTED BY USING THE JACCHIA 71/LINEBERRY 
DENSITY MODEL. 

(2) p IS THE DENSITY AT THE "GEOMETRIC" PERIGEE. 

(3) p p AND h s ARE RECOMPUTED AT EACH STROBOSCOPIC STEP. 



t 
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VERIFICATION AND CHECKOUT OF DSTROB 
(a) PSU vs. DSTROB 
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EVOLUTION OF PERIGEE ALTITUDE 6 
EFFECT OF LAUNCH DATE 
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EVOLUTION OF PERIGEE ALTITUDE GOES-E/SUSS-A 
EFFECT OF LAUNCH DATE 



H 

oo 

0) O) 
P H 
cd 

Q J3 
o 

x: u 
O ed 


cd co 

<N 


E 

X 

o 

o 

o 

in 

in 

o 

o 

CO 

00 

in 

00 

H 

H 

00 

H 


II 

II 

II 

II 

O 

CX 

^0 

0 

c 

x: 

0* 

3 



d) oo 
P O) 


x: a 

O 3 
C >3 

n 

cd Cd 
<N 


(w») aanimv 339id3d 


TIME (YEARS) 


STS-6 SBS/SUSS-D 
(a) perigee time history 

h = 296 km 



TIME (YEARS) 



STS-6 SBS/SUSS-D 
(b) perigee time history 

h = 185 km 
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NECESSARY CONDITIONS FOR IMMEDIATE 
DECREASE IN PERIGEE ALTITUDE 


• PERIGEE ALTITUDE WILL DECREASE WHEN ECCENTRICITY 
INCREASES 


• FOR AN APPROXIMATE FORMULA, NEGLECT THE FOLLOWING 

effects: 

• AIR DRAG 

• SIN I, SIN € 

• ORBITAL MOTION OF THE MOON AND SATELLITE 


• THE CHANGE IN ECCENTRICITY IS APPROXIMATELY 

a *~ rio <£ . . r. /-. \1 


7^ = MEAN LONGITUDE OF THE SUN 

-fL .36 deg/day 


- Ce>5[2(2&-jn.-i* r ) 0 ']\ 

F flF TWF SUM / 


^ .62 deg/day 


+ THIS FORMULA WAS SUPPLIED BY 
PROFESSOR JOHN V. BREAKWELL 
DEPARTMENT OF AERONAUTICS AND ASTRONAUTICS 
STANFORD UNIVERSITY 


NECESSARY CONDITIONS FOR IMMEDIATE 
DECREASE IN PERIGEE ALTITUDE (concl.) 


USING THE PREVIOUS FORMULA, IT CAN BE SHOWN THAT 
THE GREATEST INITIAL IN CREASE IN OCCURS WHEN 


(A.-_n ) 0 = - Ti 

- -rrl - - + % 


THIS RESULT IS ILLUSTRATED IN THE NEXT SLIDE, 
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Launch Window STS-8 GOES-E/SUSS-A h - 185 km 
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Launch Window STS-7 TDRS/IUS h = 185 km 
















Time of Lift-off (Hours, GMT) 

GOES-E Launch Window h D0 =185 km 



Time of Lift-off (Hours, GMT) 

TORS Launch Window h nn = 185 km 





CONCLUSIONS AND RECOMMENDATIONS 


• THE MOST RAPID ORBIT DECAY OCCURS WHEN THE SUN'S 
RIGHT ASCENSION IS 90° AHEAD OR BEHIND THE 
ASCENDING NODE OF THE TRANSFER ORBIT. 

• THE TRANSFER ORBIT SHOULD BE DESIGNED SO THAT THE 
INITIAL PERIGEE ALTITUDE IS LESS THAN 232 KM (125NM) . 

• ALLOWS ADEQUATE STAY-TIME IN TRANSFER 
ELLIPSE. 

• ALLOWS POSSIBILITY OF FUTURE ORBIT DECAY. 

• TRANSFER ORBIT DECAY LAUNCH WINDOWS CAN BE GENERATED 
FOR EACH MISSION AND INTEGRATED INTO THE. COMPOSITE 
LAUNCH WINDOW. 


SUMMARY OF THE WORK 


• A SEMI -ANALYTICAL TRAJECTORY GENERATION METHOD 
(DSTROB) HAS BEEN DEVELOPED AND TESTED. IT IS 
APPLICABLE TO GEOSYNCHRONOUS TRANSFER ORBITS. 

• ORBIT DECAY PREDICTIONS WERE COMPARED TO AIR 
FORCE TRACKING DATA. GOOD AGREEMENT WAS 
OBTAINED. 

• A STUDY WAS MADE ON THE ORBIT DECAY DYNAMICS OF 
GEOSYNCHRONOUS TRANSFER ORBITS. 

• UPPER STAGE ORBIT DECAY WAS STUDIED FOR CERTAIN 
PLANNED STS MISSIONS. 


SUMMARY OF THE WORK (concl.) 


• A PROTOTYPE COMPUTER PROGRAM WAS DEVELOPED TO 
PROVIDE ORBIT DECAY PARAMETER SCANS AND LAUNCH 
WINDOW DATA. 

• A TECHNIQUE WAS DEVELOPED FOR DETERMINING ORBIT 
DECAY LAUNCH WINDOWS FOR STS UPPER STAGE MISSIONS. 

• UPPER STAGE ORBIT DECAY LAUNCH WINDOWS WERE 
GENERATED FOR CERTAIN STS MISSIONS. 


APPENDIX B 
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A STUDY OF THE LIFETIMES OF GEOSYNCHRONOUS TRANSFER ORBITS 
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A STUDY OF THE LIFETIMES OF . 

GEOSYNCHRONOUS TRANS FEE ORBITS* 

Otis F. Graf;, Jr. aiul Alan C. Mueller^ 

This study concerns the orbit lifetimes of spent upper 
stages that are used to boost NASA Space Transporta- 
tion System payloads into geosynchronous transfer or- 
bits. A semi-analyl Leal method of trajectory compu- 
tation was developed and applied to these particular 
orbits. This method allowed very fast computation of 
orbiL lifetimes. it was found tliat gravitational 
perturbations of the sun, moon and oblate earth cause 
large changes in the perigee al.Li.Lude. Lifetimes 
may vary from a few months to many years. The follow- 
ing parameters have a strong influence on orbit life- 
times: Time of year of launch; Inertial orientation 

of orbital plane; inclination of orbital plane. The 
transfer orbits can be designed to decay within one 
year if tlte initial perigee altitude is less than 
231 km. However, titer e are restrictions on orbit 
plane placement and time of year of launch. The 
time of year restriction does not apply if initial 
perigee is lowered to 185 km, 

INTRODUCTION 

Statement of the Problem 


The NASA Space Transportation System (STS) offers new economies for plac- 
ing payloads into earth orbit. There will be new opportunities to orga- 
nisations for operating their devices in space. Thus, there will be a 
much wider utilisation of space as a resource for applications. 

At the same time, however, the Increasing number of satellites will in- 
crease the possibilities of collision in space,. In the past, most 


* Presented at the AAS/AIAA As f roily mini ios Specialist Conference, 
Provincetown, Massachusetts, June 25-27, 1979. 

' Analytical and Computational Mathematics, Inc. Houston, Texas. 
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satellites were in near earth orbit, where air draft eventually forced 
them to fall back to earth. Now, however, many satellites will be pul 
into geosynchronous orbit, using Lho. lnLori.ui Upper Stage (1US) or the 
Spinning Solid Upper Stage (SSUS). 

l'or each satellite launched into geosynchronous orbit, one or more upper 
stages will be left in the transfer orbit. These orbits have a low peri- 
gee and will intersect the orbit of the Shuttle. Also, these orbits can 
have a very long lifetime. During the ten year period 1980-1990, there 
are 164 planned geosynchronous missions using thu STS Space Shuttle 
(Ref. 1). 


The purpose of this study was to investigate thu orbit lifetimes of the 
spent upper stages that are usod to boost STS payloads into geosynchronous 
transfer orbits. The upper stage vehicles remain in the transfer orbit. 


it is ussumud that the transfer orbit has an initial perigee altiLude of 
300 km and an apogee altitude of about 35,785 km. The scale of the orbit 
Is shown in Fig,. 1. It can bo soon that the perigee altitude has a rela- 
tively small scale. Thus, a small change' In the sine or shape of the 
orbit can drastically raise or lowur perigee. 



Fig. 1 Scale of Geosynchronous Transfer Orbit 
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Tina orbit geomutry is shown in l>’ig. 2. I)u lino the following symbols: 

£ Angle between the ecliptic plane and the equatorial plane 
(23.3°). 

T Direction of the vernal equinox. 

Q Angle of the ascending node of the satellite orbit on the equa- 
tor, measured from the equinox. 

1 Inclination of the geosynchronous transfer orbit relative to the 
equator. 

a u Righ i ascension of the sun. 
b 



fig. 2 Orbit Oeometry 

Since the geosynchronous orbits are placed at aero inclination, it is as- 
sumed in this study that the transfer orbit inclination will be less than 
the Shuttle parking orbit inclination. That is, I will be equal to or 
less than abouL 28 degrees. 

Summary of Results 

It was found that the gravitational perturbations of the sun, moon and 
oblate earth can make large changes in the, perigee altitude of the geo- 
synchronous transfer orbits. Some, transfer orbits may have u lifetime 
of only a few mouths, in uLher cases, the perigee altitude is raised so 
that the orbit will have an extremely long lifetime. 

'L'he important results are listed below: 

1. The following parameters have a strong influence on orbit lifetimes: 
Time of year of launch. 

Inertial orientation of orbital plane. 

Inclination of orbital plane. 
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2. The moat severe decay in perigee altitude occurs when the initial 
longitude of ascending node is near 180 degrees. 

3. If the initial perigee alLiLudu is less than 231 km (125NM), the 
transfer orbit can be designed to decay within one year. 

4. Geosynchronous transfer trajectories may have lifetimes of many 
years for some cases. 


EFFECTS OF THE DIFFERENT ORBIT PERTURBATIONS 


The geosynchronous transfer orbits are strongly perturbed by the follow- 
ing forces: Gravitation of the oblate earth (J 9 ), Gravitation of the 

moon (point mass), Gravitation of the sun (point mass), Atmospheric drag. 

All four of the above perturbations are of equal importance. There is at 
present, no analytical solution Lo the motion of a satellite in a geosyn- 
chronous transfer orbit (Refs . 2 and 3). Since all perturbations are of 
about the same magnitude, a linear perturbaL ion method does not yield an 
accurate solution. Therefore, trajectory generation must he done with 
numerical or seml-analyl ieal methods. This aspect of the problem Is dis- 
cussed in the next section. 

The effects of the different gravitational perturbations on perigee alti- 
tude are shown in Fig. 3. The small, oscillations are caused by the 
monthly revolution of the moon in its orbit, it can be seen- in Fig. 3 
that the magnitudes of the perturbations are all about the same. 
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It was found that atmospher i c drag duos not have an important of feet on 
orbit decay when perigee ailiiude is above IbO km. The effects of air 
drag on orbit decay are shown i n Kip,. ^ . The drag model was based on 
the burn-out weight and average cross-sectional area of the Spinning 
Solid Upper Stage (SSUS). As the gravitational perturbations bring peri- 
gee down, air drag causes the apogee to drop. As the apogee comes down, 
the sun and moon perturbations on perigee are less severe. Therefore, 
perigee drops less rapidly Lor Lhe drag-perturbed case than for the case 
whore drag is neglected. 



Fig. 4 1CI loots ol Drag on Orbit Decay 


For this study it was important to determine r Lhvcsliold pel'igee attitude } 
defined as the maximum perigee altitude for which the satellite will deor- 
bit within 90 days. That is, at any perigee altitude below the threshold, 
atmospheric drag will be certain to deorbit the satellite within 90 days. 


Tiie determination of the threshold perigee alLitude 


h , was done by 
pt 


taking a case where the gravitational perturbations had the strongest 

effect on raising perigee. .Several trajectories were simulated using the 

KSFAST program on the UNiVAC 1110 computer (Ref. 4). Each trajectory was 

initialized at a different perigee altitude h . The evolution of 

po 

perigee altitude h Cor several cases is shown in Fig. 5. It was 


at altitudes 


found that orbit lifetime was extremely sensitive to h 
near 100 km. For h < 101.5, the orbit lifetime was less than 90 days. 
Therefore, h was chosen to be 100 km. This value is in agreement 
with the remarks made in Section 5.1 of Ref. 3. 



Fig. 5 Effects of Drag on Low Perigee Orbits 

In the remaining discussion, it is assumed that when the perigee altitude 
drops to 100 Ion, the orbit will rapidly decay. This threshold altitude 
will be important in the discussion of orbit lifetime parameter scans. 

A SEMI-ANALYTICAL TRAJECTORY GENERATION METHOD 

Trade-Off Studies with Existing Routines 


In order to get a global picture of the stability of Lhu geosynchronous 
transfer clqss of orbits, their lifetimes must be determined as a func- 
tion of the problem variables. This will require that many hundreds or 
thousands of trajectories be generated on a set of initialisation para- 
meters. These orbit lifetime parameter scans will give a general picture 
of how orbit decay is affected by a variety of control variables, such as 
time of launch, initial orbit placement, etc. 
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The complete soL of L raj eel ory programs a v;i i lab I e at the Johnson Space 
Center (JSC), Mission Planning and Analysis Division (MPAD) , were inves- 
tigated for their application to this prohJem. existing analytical pro- 
grams do not have an accurate representation of the perturbing forces. 
Therefore, numerical step-by-step programs were investigated. The required 
execution times (per revolution) of most numerical programs eliminated 
their application to this problem. Two MPAD programs found to be nearly 
applicable were KSFAST (Ref. 4) and STEPR (Ref. 5).' Both programs are 
available on the UN1VAC 1110 computer system. 


The execution times of STEPR and KSFAST are shown in Table 1. 
case was a two year trajectory generation of a geosynchronous 
satellite. STEPR was found co be much faster than KSFAST. 


The test 
transfer 


Table I 


COMPUTER EXECUTION TIMES FUR A TWO YEAR TRAJECTORY 


STEPR 

100 see. 

UNI VAC 

1110 

KSFAST 

HU0 sec . 

UN L VAC 

1110 

DSTROB 

1 . 5 sec . 

UN1VAC 

1108 


A typical scan over the 360 degrees range of an angular variable would 
require 36 trajectories, if the scan interval, was 10 degrees. The ini- 
tial longitude of the ascending node S'2 is a typical angular variable 

for this problem. Computer execution times for such a scan are shown in 
Table 2. Each trajectory in the scan was propagated lor two years. 

Table 2 


COMPUTER 

EXECUTION 

TIMES FOR A 360 DECREE 
(36 points) 

SCAN 

STEPR 

60 

min. 

UNIVAC 

11.10 

KSFAST 

480 

m in . 

UNIVAC 

1110 

DSTROB 


.9 min. 

UNIVAC 

1108 


The data in Table 2 demonstrates that both STEPR and KSFAST would be un- 
suitable for parameter scans. Therefore, no available JSC program could 
be used fur an important part ol this sludy. Fur ibis reason , ACM imple- 
mented and used the semi-analytical algorithm described below. 

T he DSTROB Algorithm 

The semi-analytical algorithm in DSTROB is based on the stroboscopic 
method of orbit generation. This method has been extensively developed 
by Roth (Refs. 6, 7 and 8). The stroboscopic, method is particularly 
suited to high eccentricity orbits. Roth has applied the method to 
launch window studies for orbits with eccent ric Ltes near 0.95 (Ref. 7). 
Janin discusses the decay of such orbits (Ref. 9). 

(7) 


aaaaaaaaaaiatvt . ■ ■ 


1 


The DSTROB algorithm makes use of Delaunay-Similar (DS) elements (Ref, 10) 
and the first order stroboscopic method described in Refs. 6 and 7. 

(Hence the name DS-STROB, or DSTROB) . The perturbing forces include: 
Earth's oblateness (Jg) , Sun point mass, Moon point mass, and Atmospheric 

drag. The sun and moon gravity models include only the first term in the 
distance ratio expansion. Also, the sun and moon are held fixed over one 
revolution of the satellite. The drag model is taken from King-Ilele 
(Ref. 11) where the scale height and the density at perigee are found 
from a atmosphere model which includes the solar and seasonal variations 
in density (Ref. 12). 

Comparisons between DSTROB and two numerical methods are shown in 
Rigs. 6(a) and 6(b). Case 1 (Rig. 6(u)) is a comparison with a method 
that integrates numerically the Poincare-Similar elements (PSu). For 
Case 2 (Rig. 6(b)), the perigee curves l'oi* KSRAST and DSTROll are so simi- 
lar that they appear as one curve. The deviation in perigee altitudes in 
Case 1 occurs because DSTROB does not account for the long period or 
second order perturbations in the eccentric i Ly due to the zonal harmonics. 
This slight deviation in perigee over an extended period of time re, suits 
in Lhe differences in h predicted by Lite Lwo methods. 



Rig. 6 Verifieac Lon and Check-out of DSTROB 
(a) PSu vs DSTROB 
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Fig. 6 Verification and Check-out of DSTROB 
(b) KSFAST vs DSTROB 


Tests such as those shown in Fig. 6 indicate that DSTRUB gives good ac- 
curacy for orbit predictions of about two years. Computer execution 
times for DSTROB are shown in Tables 1 and 2. It j.s seen that DSTROB is 
much faster than the other two programs. An additional feature of DSTROB 
is that it can be executed in the demand (interactive) mode on the 
NASA/ JSC UNIVAC 1108 Computer system. This greatly facilitated the data 
generation and analysis part of the study. 

EVOLUTION OF PERIGEE ALTITUDE 


This section will illustrate some cases of the long term variation in 
perigee altitude. As discussed in the previous sections, the gravita- 
tional perturbations can have a significant effect on perigee altitude. 
Thus, perigee can be either raised or lowered, depending on the orbital 
initial conditions. 


Examples of the evolution in perigee altiLude over two years are shown in 
Figs. 7 and 8. These cases were taken from the COES-E mission that is 
planned to be launched from Space Transportation System (STS) Flight 8. 
The upper stage booster will be a SUSS-A in the following transfer orbit'. 
Apogee Altitude (h ) 38,034 Ion 


Perigee Altitude (h ) 

P 

Right Ascension of Ascending Node (Q) 
Argument of Perigee (u) 

Launch Date 


296 km 


350° 

180° 

26 March 1981 
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For the drag model, It is assumed that the SUSS-A burn-ouc weight is 
316 kg, with an average (tumbling) cross-sectional area of 3,92 square 
meters. The coefficient of drag is assumed to be 2.2. 

Two different launch dates are considered in Figs. 7. In Figs. 8, differ- 
ent values of the initial ascending node are considered. For each of 
these comparisons, three initial perigee altitudes are investigated: 

296 km, 231 km and 185 km. The different perigee altitudes were included 
since NASA is giving some consideration to launching the satellites into 
transfer orbits with lower perigee altitudes than currently planned. 



Fig. 7 Evolution of Perigee Altitude 
(Effect of Launch Date) 

(a) h =296 km 
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Fig. 7 Evolution of Perigee AitiLude 
(Effect of Launch Date) 

(c) h = 185 km 
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Fig. 8 Evolution of Perigee Altitude 

(Effect of Orbit Plane Placement) 

(a) = 29b km- 


/~Q 0 = 1^0“ 

200 f \ 

\ ^ ' 

•N J J 

\ Vv-/ 

160 \ /J 


I II I I I il I L'llllll I | I HUM 


Ii - ,'iH,0;M km 


■Q a = 250° 


I. mint'll Diitu. 20 Miil'cl) 1981 
V'~\a 


^ e, v 

V' l . a 


TIME (YEARS) 

Fig. 8 Evolution of Perigee Altitude 

(Effect of Orbit Plane Placement) 
(b) h =231 km 
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Fig. 8 Evolution of Perigee Altitude 

(Effect of Orbit Plane Placement) 
(c) h = 185 km 
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It is seen that launch date and ascending node have a significant effect 
on evolution of h , and hence orbit lifetime. (This is further illus- 


trated by the parameter scans in a later section.) However, with the 
nominal initial perigee altitude (h = 296 km) none of the orbits in 

po 

this study were found to decay within two years. Only by going to lower 
altitudes (hp Q <, 231 km) , did the transfer orbits decay for certain ini- 
tial conditions. 


The following are some additional observations on Figs. 7 and 8: 


1, As the apogee comes down, the perturbation of the moan is less severe 
and the small oscillations begin to smooth out. (Figs. 8(b) and 8(c)) 


2 * 


With hp Q = 185 km, there are several launch window opportunities such 
that the transfer orbit decays within a few months. 


3. The DSTROB routine did not simulate the satellite's destruction, but 
only carried the trajectory to the point where reentry would occur 
within a few revolutions. 


A. For the planned launch conditions of CEOS-E/SUSS-A, the transfer or- 
bit will not decay within two years, even if perigee is initialized 
at 185 km (Figs, 7(a), 7(b) and 7(c)). • 
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EXAMPLES OF ORBIT DECAY FROM AIR FORCE TRACKING DATA 


The previous section presented computer simulations that depicted large 
fluctuations in perigee altitude. The question arises as to whether 
these simulations agree with the observed decay of satellites in geosyn- 
chronous transfer orbit. In this section the computer simulation (the 
DSTROB program) is compared to satellites that were tracked by the U.S. 
Air Force. 

Trajectory data on all observed satellites are contained in a periodical 
called CLASSY (Ref. 13), published by the North American Air Defense 
Command (ADCOM) . The satellite orbit data presented in this section were 
taken from these documents. Each of the four satellites presented here 
had orbits that were similar in apogee, eccentricity and period to that 
of a geosynchronous transfer orbit. 


Comparisons of DSTROB with Air Force tracking data are shown in 
Figs. 9(a) - 9(d). The solid line in each figure is the output from 
DSTROB. The triangles represent data that were taken from precise or- 
bila] elements, obtained via private eummun ieal ions with ADCOM. These 
precise elements were used as inputs to DSTROB. When two triangles ap- 
pear in a figure, only one was used to obtain Lhe input vector to DSTROB. 



Fig. 9 Comparison of DSTRO 1 ' with Air Force Trucking Data 
(a) Satellite No. 5713, J = 05.1° 
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Fig, 9 Comparison ol DSTRUI1 With Air Force Track Lug Data 
(d) Satellite No. 8?v) r \ 1 ui,l) u 

Vi 

Thu dots in Figs. 9(a) - 9fd.l represent Jala taken Item the orbital e lu- 
men t a published in the CLASSY document a. The epoch lor a act of CLASSY 
elements is known to be within a three monLh interval of time. Wo at- 
tempted to refine the epoch by starting, from the precise elements and 
using an average value of the mean motion. Therefore, the time coordi- 
nates for the dots can have an error ol 1-) months. Drag model data was 
not available on the four satellites, There l ore, atmospheric drag effects 
could not be simulated. Even so, the DSTEUlt program gives quite good 
agreement with the Air Force data. Satellites No. f> 7 1 3 and No. 8792, decay 
rapidly duo to gravitational perturbations, '(’heir orbits are unstable be- 
cause of the large inclinations. Fig. 9(d) is another example where the 
purigoe altitude decays less rapidly (for awhile) when drag becomes sig- 
nificant. 

ORlirr LIFETIME PARAMETER SCANS 

The remaining part of this study is to del ermine orbit l Hotline as a 
function of the, free, parameters of the problem. The parameters of inter- 
est are the set of initial orbital elements, defined as follows: 
a - semi-major axis, 

e * eccentricity 

1 = inclination 

U) » argument of perigee, 

Q a longitude of the ascending node (Fig. 2) , 

T - epoch. 

Oh) 


nil Lai values lor 4 and e arc spoiitieU in the mission objectives. 

It w J1 also he cogui. red that pci tgee Is ihmi | hr oguatoi and lli.il t | u . 
lncUiuuum be near JiH.ft , Thoroloro, (j .uul I aro mil soluble scan 


pa rameters, 


lho remaining parameters oi interest are: 

initial longitude of Lho ascend ing node. 
a date oL launch, 

lho Lime of day of launch can ho related to Q , There! ore, Tj aliould 

sunnily Lho day ol year ot launoli. A physical parameter nilgliL bo mot‘o 
usofui. 1 lhuo, Lho right asoemslon ol lho huh a Iwiv Kip,. 2) was used 
tor Lho time ol year pa fame tor, ° 

In piaoliee^ 1! Ifi geosynoln ouous mission.-. involve a vaiiety ol launoli eon- 
st minis. Those oonsL ra i n t s depend, in pail, on the mi mi Lon and payload. 
XL was, there tore, oulsi.de the scope ol this study to include, miss ion 
cun si mints on geosynehrouous t i ansi or oihit hlelimo. 

The pa ramo Lor scans wore done by oompul ing in one year trajectories ol a 
geosynchronous transler satellite. It was assumed that the launch occur- 
red during 191)0. For each trajectory, lho minimum perigee in one year 
was determined, The angular parameter was incremented by 10 degrees over 
the JtiO degree range. 
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Scans on the Ascend ing Nod e 


The parameter of interest, in Litis case was Lite initial iungilude oi Lite 

ascending node O , which deseribes Lhc inertial placement of the or- 
o 

bital piano relative to Lite pianos of Lite sun, moon and equator. This 
paramo tar has an important of fee l on orbit evolution. 

Plots of minimum perigee as a function of Q are. shown in Figs. 10(a) - 

10(c). Three different initial perigee altitudes were, investigated: 
hp Q s 185 km, 231 km, 278 km. For each ease, trajectories were initial- 
ized for summer (June), winter (December), spring (March) and fail 
(September). The epoch was 1980. 

The lowest perigee altitudes always ooeur near Q = 200°. Trajectories 
that are initialized in summer and winter allow lower perigee altitudes. 

Scans on the Launch Da te 

The time, of year paramter was taken to he I he dillereueo between the 
sun's right ascension at launch and Lite initial Longitude of the ascend- 
ing node, i.e., Ct^-O < both angles are defined in Fig. 2. This 

parameter was chosen because if describes the position of the sun in its 
orbit relative Lo the satellite orbital pi. me. 



<*. " Q q lutu.) 

11 Orbit Lifetime Parameter Scans 
Time of Year ul Launch 


(a) h p£) = 185 km, 


l u - 28 . if 


£> = 180 L 


Parameter scans are shown in Figs. 11(a) and 1.1(b) lor an initial perigee 
altitude of 185 km, for the two eases = 180° and Q - 0. As would 
be expected from the results of Pig. 10, Llie case Q ^ = 180 gives a con- 
sistently low minimum perigee, regardless of the time of year. 

For other values of Q q , the perigee is lowest when the sun begins 90° 
(or 270 ) ahead of the ascending node. 



Fig. II Orbit Lifetime Parameter Seans 

Time of Year ul Launch 

(b) li = 185 km, 1 = 28.6°, Q = 0 
po o ** o 

Discussion of Results 

1. In Figs. 10 and 11, the threshold perigee altitude is drawn as a 
solid line at 100 km. It can be assumed that when the minimum peri- 
gee altitude falls below this line, the orbit will decay within one 
year. 

2. The parameter scans shown in this section are for I = 28.6°. We 

found that the results were similar for inclinations less than 28°. 

However, these figures should be used with caution if I is much 

different from 28°. We intend in the future to publish a catalog of 
parameter scans for various epochs and 'inclinations. 
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If h - 185 km, there will bo Launch opportunil Los throughout Llu> 

year such that the transfer orbit wiJ I decay within one year 
(Pigs. 10(a) and 11 (a)). For h - -31 km, good chances Cor orbit 

decay occur only Cor summer and winter launches (Fig. 10(b)). 

If h - 278 Ion (150 N.M.) there is little possibility of rapid de- 

po • 

cay of the transfer orbit. 
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